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Abstract 
^^ , We extended the study of relativistic effects in neutron-deuteron scattering to the exclusive breakup. To 

C^ . this aim we solved the three-nucleon Faddeev equation including such relativistic features as relativistic 

kinematics and boost effects at incoming neutron lab. energies _E^"'' — 65 MeV, f56 MeV and 200 MeV. 
As dynamical input a relativistic nucleon-nucleon interaction exactly on-shell equivalent to the CD Bonn 
potential has been used. We found that the magnitude of relativistic effects increases with the incoming 
neutron energy and, depending on the phase-space region, relativity can increase as well as decrease the 
nonrelativistic breakup cross section. In some regions of the breakup phase-space dynamical boost effects 
(«r^ ' are important. For a number of measured exclusive cross sections relativity seems to improve the description 

o 

O . PACS numbers: 21.45. -Hv, 24.70.4-8, 25.10.+S, 25.40.Lw 

O 






(N 

^ . 
Oh. 

< 



of data. 



X 



I. INTRODUCTION 

The high precision nucleon-nucleon (NN) potentials which describe very weh the NN data set 
up to about 350 MeV [l|, y, y| form a very firm basis for a study of three-nucleon (3N) reactions. 
Powerful computers and development of modern algorithms provided numerically exact solutions 
of the 3N Schrodinger equation both in the momentum and coordinate representation. This per- 
mitted theoretical calculations for cross sections and spin observables in elastic nucleon-deuteron 
(Nd) scattering and breakup processes with different dynamical assumptions about underlying 
nuclear forces 4]. With increasing amount of precise 3N elastic scattering data it turned out, 
that nonrelativistic description based on pairwise forces only is insufficient to explain the data at 
higher energies of the 3N system. Adding a 3N force (3NF) to the pairwise interactions led in 
many cases to a better description of the data |5|, la, M> M> M> UM> \lM UA, \1^ \1M> ll^l- However, at 
energies higher than ~ 100 MeV current 3NFs only partially improved the description of data, 
leaving in some cases discrepancies which were comparable in magnitude to the effects of the 3NFs 
themselves |d, li Id, Ie Ld , lul, ui ImuAuM 



This situation triggered investigations of the 3N continuum taking relativistic effects into ac- 
count. In ref. jlQ] we performed such a study for Nd elastic scattering. We extended the Hamilto- 
nian scheme in equal time formulation worked out in 17] to the 3N scattering taking as a starting 
point the Lorentz boosted NN potential which generates the NN t-matrix in a moving frame via a 
standard Lippmann-Schwinger equation. The NN potential in an arbitrary moving frame is based 
on the interaction in the two-nucleon cm. system, which enters a relativistic NN Schrodinger or 
Lippmann-Schwinger equation. The relativistic equation differs from the nonrelativistic one just 
by the relativistic form of the kinetic energy. We constructed the relativistic two nucleon (2N) po- 
tential by performing an analytical scale transformation of momenta, which relates NN potentials 
in the nonrelativistic and relativistic Schrodinger equations in such a way, that exactly the same 
NN phase shifts are obtained by both equations |18l |. 

In our first study 16] we looked for changes in elastic neutron-deuteron (nd) scattering observ- 
ables when the nonrelativistic form of the kinetic energy was replaced by the relativistic one and 
a proper treatment of boost effects and Wigner rotations of spin states was included. It turned 
out, that the effects of spin rotations in the studied energy range up to 250 MeV were practically 
negligible for elastic scattering cross sections and analyzing powers. The relativistic effects for the 
elastic scattering cross section were significant only at higher energies and restricted to the very 
backward angles, where relativity increased the nonrelativistic cross section. The decisive role was 



played by the boost effects which reduced the transition matrix elements at higher energies and 
led, in spite of the increased elastic scattering relativistic phase-space factor as compared to the 
nonrelativistic one, to rather small effects in the cross section. 

In view of upcoming measurements of the higher energy Nd breakup process we would like to 
extend the study of ref. [la] to this reaction. The first results presented in [ig| revealed a unique 
selectivity of the complete breakup reaction useful to investigate the pattern of relativistic effects 
and to study their significance. Here we would like to study how relativistic effects are distributed 
over the breakup phase-space and if and to what extent the existing data verify the predicted 
effects. 

The paper is organized as follows. In Sec. II for convenience of the reader we shortly explain 
the relativistic features underlying our treatment of a relativized Faddeev equation in the 3N 
continuum. The very detailed presentation which incorporated the definition of the boosted two- 
body force, the various two-and three-body states in general frames, the WijEjner rotations and the 
singularity structure of the relativistic free 3N propagator were given in |1^]. In Sec. Ill we give 
relativistic and nonrelativistic formulas for the transition matrix elements and the breakup cross 
section. In Sec. IV we apply our formulation using a relativistic NN interaction which is exactly 
on-shell equivalent to the nonrelativistic CD Bonn potential and solve the relativized 3N Faddeev 
equation with different approximations for the boost. At two incoming neutron energies, 65 and 
200 MeV, we look for the magnitude and the distribution pattern of relativistic effects over the 
entire phase-space of the breakup reaction. For the cases where the data are available, we compare 
them to theoretical predictions. Finally, Sec. V contains a summary and outlook. 

II. FORMULATION 

When nucleons interact through a NN potential V, the breakup operator T satisfies the Faddeev- 
type integral equation 4, l2i| 

T\^ > = tP\^ > +tPGoT\^ > . (1) 

The 2N t-matrix t results from the interaction V through the Lippmann-Schwinger equation and 
the permutation operator P = P12P23 + ^13^23 is given in terms of a transposition Pij which 
interchanges nucleons i and j. The incoming state \(j) >= [^0 > \4'd > describes the free nucleon- 
deuteron motion with relative momentum % and the deuteron wave function \(j)d. >■ Finally Go is 
the free 3N propagator. 



The full nd breakup transition operator Uq is given in terms of T by 



2C| 



2l| 



Uo = (i + P)r. (2) 

This is our standard nonrelativistic formulation, which is equivalent to the nonrelativistic 3N 
Schrodinger equation and respects the boundary conditions. The formal structure of these equa 
tions in the relativistic case remains the same but the ingredients change. As explained in 
the relativistic 3N Hamiltonian has the same form as the nonrelativistic one, only the momentum 
dependence of the kinetic energy changes and the relation of the pair interactions to the ones in 
their corresponding cm. frames changes, too. Consequently all the formal steps leading to Eqs.(^ 
and remain the same. 

The relativistic kinetic energy of three equal mass (m) nucleons in their cm. system can be 



written in terms of the momentum dependent 2N mass operator 2u){k) = 2\J m? + k"^ and the 
momentum of the third nucleon q as 



Ho = ^{2oj{k)Y+q^ + ^m? + q^. (3) 

Here k and —k are the momenta of two nucleons in one of the two-body cm. subsystems and —q 
is the total momentum of this chosen two-body subsystem. Any of the three possible two-body 
subsystems can be taken. 

The boosted 2N potential in the 2N frame moving with momentum q is taken as J22] 



V{q) = ^{2uj{k)+vY + q^-^{2i^{k)f+q^, (4) 

where V{q) for g = reduces to the potential v defined in the 2N cm. system. Note that also 
in that system the relativistic kinetic energy of the two nucleons has to be chosen, which together 
with V defines the interacting two- nucleon mass operator occurring in Eq.(|3J). 

The boosted 2N t-matrix t{k, k ']q) fulfills the relativistic 2N Lippmann Schwinger equation, 
which in a general frame reads 



'{2u;{k') 2 + g2_^(2tj(A:") ■^ + q^ + ie 

Using Eq. Q) the relativistic 2N Schrodinger equation for the deuteron in a moving frame can 
be cast into the form 

Mk) = , ^, ^ fd'k'Vik,k';qo)Mk'), (6) 

'Mj + q,' - J{2u:{k)Y + q,^ ^ 



where Md is the deuteron rest mass and \/ Mj + ^q^ its energy in the moving frame. 

The new relativistic ingredients in Eq.(^ will therefore be the boosted t-operator and the 
relativistic 3N propagator 

where Hq is given in Eq. (3) and E is the total 3N cm. energy expressed in terms of the initial 
neutron momentum rfn relative to the deuteron 



E = ^{MaY + q^^ + ^m^ + q^\ (8) 

We solve numerically Eq.(^ in its nonrelativistic or relativistic form for any NN interaction 
using a momentum space partial wave decomposition. In the nonrelativistic case the partial wave 
projected momentum space basis is taken as \pqa >= \pq{ls)j{X^)IJ{t^)T >, with the magnitudes 
p and q of standard Jacobi momenta (see ^j) and {ls)j two-body quantum numbers with obvious 
meaning. The quantum numbers (A1/2)I refer to the third nucleon (its motion is described by the 
momentum q), and J is the total 3N angular momentum. The subsystem isospin t couples with 
the spectator isospin to the total 3N isospin T. 

In the relativistic case the nonrelativistic relative two-nucleon momentum p is replaced by k, 
where 2N cm. momenta k and —k are related to general momenta of these two nucleons, say p2 
and p3, by a Lorentz boost: 

k = k{p2,P3) 

/ 2 — E^ 

2 E2 + E3 + J{E2 + E3)^-pi 



with Ei = Jm? + p^"^ and p23 = P2 +P3- The third spectator nucleon has momentum pi which 
together with the total two-nucleon momentum p23 adds up to zero in the 3N cm. system. It is 
the momentum q = pi which replaces the nonrelativistic Jacobi momentum q in the relativistic 
case to describe unambiguously a 3N configuration. 

The construction of the momentum space partial wave basis in the relativistic case starts 
with definition of the 2N subsystem partial wave state \{ls)jkfi;0 > defined in the 2N cm. 
subsystem [l6|. Boosting this state to the 3N cm. system along the total momentum p23 
and coupling it there with the state of the spectator nucleon leads to a 3N partial wave state 
\k q = pi a >= \kpi{ls)j{\\)I{jI)JM > \{t}2)TMT >■ Since the performed boost generally is 
not parallel to the momenta k and —k of the two nucleons in their 2N cm. system it leads to 



Wigner rotations of their spin states. This rotation comphcates the evaluation of the partial wave 



m 



it was shown 



representation of the permutation operator P in the basis \k q = pi a >. In 

that this representation can be given in a form which resembles closely the one appearing in the 

nonrelativistic case jj, l23| | 

6{ k-TTi) 6{k' - TTa) 

Gaa'iq,q,x), (10) 



< k q a\ P \k' q' a' > = / dx 



Ni{q,q',x) N2{q,q',x 
what leads after projecting Eq. (^ onto \kqa > to an infinite system of coupled integral equations 



analogous to the nonrelativistic one 



20|: 



PCX) PL 

< kqa\T{E)\(j) > = < kqa\tP\(l) > +Y.Y. dq'q'^ dx 



"1 < kQt^'^XE - Vm'' + q')\^iia > 



^ Ni{q,q',x) N2{q,q',x) tts'-' 

x^^^^ . (11) 

Xq + l€ — X 

The geometrical coefficients Gaa'iq,q' ,x), the coefficients Ni{q,q',x) and N2{q,q' ,x), and the 
momenta tti and 7r2 stem from the matrix element < kqa\P\k'q'a' > of the permutation operator 
(Eqs. (C6-C8) in ref 16]). The last part ^ :i^_x results from the free propagator (Eqs. (34-35) in 



ref. |l^). The quantum numbers in the set a differ from those in a only in the orbital angular 
momentum / of the pair. 

The main difficulty in treating Eq. (|llj) is caused by the singularities of the free propagator Go 
which occur in the region of q and q' values for which \xq\ < 1. In addition, at q = qo there is the 
singularity of the 2N t-matrix in the '^S'l — ^ Di partial wave state, where the deuteron bound state 




la,|2i|. EquationjUU is solved by 

201. 



exists. How to treat those singularities is described in detail in 

generating its Neumann series, which is then summed up by the Pade method 

Due to its short-range nature, the NN force can be considered negligible beyond a certain value 
jmax of the total angular momentum in the two nucleon subsystem. Generally with increasing 
energy jmax will also increase. For j > jmax we put the t-matrix to be zero, which yields a finite 
number of coupled channels for each total angular momentum J and total parity vr = { — Y^'^ of 
the 3N system. To achieve converged results at our energies we used all partial wave states with 
total angular momenta of the 2N subsystem up to jmax = 5 and took into account all total angular 
momenta of the 3N system up to J = 25/2. This leads to a system of up to 143 coupled integral 
equations in two continuous variables for a given J and parity. 

As dynamical input we used a relativistic interaction v, which is defined as partner of the 
relativistic kinetic energy, generated from the nonrelativistic NN potential CD Bonn according to 



the analytical prescription of ref. |l2]- In ref. 1T\ it was shown that the explicit calculation of the 
matrix elements V{k, k '\q) according to Eq. (4) for the boosted potential requires the knowledge 
of the NN bound state wave function and the half-shell NN t-matrices in the 2N cm. system. 
Here, as in our study of the elastic nd scatte ring |16l |. we do not treat the boosted potential matrix 
element in all its complexity as given in ref. [13] but restrict ourselves to the leading order term in 
a q/io and v/lo expansion 

V{k,k'-q) = v{k,k')^[l - ]■ (12) 



S\lm? + k "^sjim? + {k'Y 

In order to study importance of the boost effect we will present in addition to this approximation 
also the results for two more drastic approximations. In the first one the boost effects are neglected 
completely 

V{k,k';q) = v{k,k'), (13) 

and in the second one the k-dependence of the first order relativistic correction term is omitted 

Vik,k';q) = v{k,k') (l-^\. (14) 

The quality of these approximations can be checked by calculating the deuteron wave function 
(j)d{k) for the deuteron moving with momentum q using Eq.®. This wave function depends 
only on the 2N cm. relative momentum k inside the deuteron and is thus independent of the 
boost momentum q. When the boost effects are fully taken into account the solution of Eq.Q 
must provide exactly the deuteron binding energy and the D-state probability equal to the values 
for the deuteron at rest. We checked in [1^1, that neglecting the boost totally or omitting the 
k-dependence of the first order term are poor approximations, especially at the higher energies 
(250 MeV) we studied. In contrast, the approximation given in Eq. (|12|) appears acceptable, even for 
the strongest boosts, reproducing closely the deuteron binding energy and the D-state probability 
for the deuteron at rest. Relying on that result we have chosen the expression (|12p for the boosted 
potential in the following investigations. Since the solution of the 3N relativized Faddeev equation 
including Wigner spin rotations due to complicated form of G^a' in Eo. HlOj) requires much more 
computer time and since we restrict in this study to the breakup unpolarized cross sections only 
we neglected the Wigner rotations completely. 



III. RELATIVISTIC AND NONRELATIVISTIC BREAKUP CROSS SECTION 

The exclusive breakup measurements d{n, Ni A^2)-^3 are performed in the lab. system with two 
of the three outgoing nucleons (A^i and A''2) detected in coincidence by detectors placed at fixed an- 
gles (O'f^, (j)^i^) and {d'f^^ (j)2''^) and their kinetic energies E^i^ and EI^'^ measured. The experimental 
events are then located in the E-^ — E2 energy plane along a kinematical curve determined by 
the energy and momentum conservation. For each point on this curve the nucleons have definitive 
momenta. The breakup observables are normally shown as a function of an arc-length S of that 
kinematical locus (the starting point of which is defined according to some convention) or as a 
function of energy E^^^ . Theoretical predictions for different observables are obtained from the 
matrix elements of the breakup transition operator (EqEJ 

< MU^\<t> > = <(po\{l + P )T\^ >, (15) 

where the state \<J)q >= \piqimim2m'iViV2V^ > describes the relative motion of free outgoing nucle- 
ons and specifies their spin (rrij) and isospin (fj) magnetic quantum numbers. This relative motion 



is described by standard Jacobi momenta {pi,qi) [J, |23| which are given in terms of individual 
momenta ki of the three nucleons in a particular kinematically complete breakup configuration by 
i,j, k = 1,2,3 and cyclic permutations 

Pi = ^(^3 -kk) 

2 ^ 1 ^ 
Qi = ^{ki-^{kj + kk)). (16) 

Applying in Ea. HlSj) the permutation operator to the left provides three contributions to the 
breakup amplitude 

I < (po\{l + P)T\(j) > = 1 < piqimiMi\T\(l) > +2 < piqimii^i\T\(p > +3 < piqimiUi\T\(j) > 

= 1 < 'Piqi'mim2m:ih'ih'2iy3\T\(t> > 

1^ 3^ ^ 1^ 
+ 1 < -^P^ - 7^1'Pi - ■^qi-,'m2m3mii/2iy3J^i\T\cl) > 

13 1 

+ 1 < --pi + -qi,-pi - -qi,m3mim2h'3h'ih'2\T\(i) >, (17) 

where the subscripts on the left side of the matrix elements indicate the leading nucleon. 

For the nonrelativistic case invariance of Jacobi momenta and amplitudes < pqa\T\(l) > permits 
calculating observables directly in the 3N lab. system. However, in the relativistic case the am- 
plitudes < kqa\T\(j) > are provided in the 3N cm. system and the transition matrix elements are 
first calculated in that system. In order to compare with the data measured in the lab. system, 
the proper Lorentz transformation from the 3N cm. to the lab. system must be performed. 



In the 3N cm. system each term in Ea. (|17(l is given by 



I < Piqi\T\(t) > = —^^-<k{k2,h),h\T\(l)>, (18) 

N{k2,k3) 



where 



^ (19) 



^■3)2 



\\{E2 + E3)^{E2 + E3)^-{k2 + k 

Assmning that the incoming beam of neutrons with the spin projection /u moves along the z-axis, 
the matrix element < k{k2, k^), ki\T\(l) > follows from the calculated amplitudes < kqa\T\(/) > by 

I s j 



<^(fc2,fc3),^i|r|0>=5:^^^'^+'"-™^-™^-™^(fc,^i)5: Vi^L^ A i / 

J,7r IXL jsItS 

[l S J 

{LSJ\fi + ma - mi - 1712 - m^-.rni + m2 + m^,^ + md){s-S\m2 + ms, mi, mi + m2 + ms) 
{t-T\u2 + U3: -^)h^t\u2, Z^3, 1^2 + i^s) < k ki a\T\(J3 >, 

where rrid is the spin projection of the deuteron. 

The 3N cm. breakup amplitude < 0ol(l + P)T\4' > provides unpolarized cross section in this 
system. Thereby the exclusive relativistic cross sections along S-curve ( ^^ ^^ dg ) "^ "-"^ projected 
on the El axis {^^^-^^)i^-^-) are given by 

and 

( ^'" )(-■) = (2vr)^ J- V 1 < 0o|(l + P)T\<^ > piS:!!f!:!L!!^, (21) 
^dnid^2dEi' ^ 2-3„^^^ ^"'^ ^ '^ ^lEf^E'^^'Es \B\ ^ ' 

with the (invariant) flux I 



T — I "1 "2 I — ^ (-pinTpin f^ in , ff in\2 (rnin,,^in\2 (00\ 

^ = \^ ~ ;^l - E^n^in V i^l -^2 -:Pl -^2 i - l™"! "^2 J ' l^^j 

and 

El P3 ■ pi 



A = 1 



B ^ I- ^^^^. (23) 



E3 pI 

E2 P3 ■ P 2 

E3 vl 



The lab. cross sections follow from the 3N cm. cross sections by the corresponding Jacobians 

/ d a y^fc _ . d a .^^ 
^dQidQ.2dEi ^dnid^2dEi 



j5 ^5 ( _dS_\c.m. 

^dnidn2dS' ^d^id^2dS' ' (g)''^''- ^ '^ 



with 



T'crn. ^c.m. 



J 



Ink I -^3'™' _ _Eg-— P3 '"- ^c.m. . -cm. 



dS v/I^Tb^ P^j 



dEi \B\ 

In Eqs. (|2UI25p m™' are the rest masses of the incoming particles and E"^" (Ei) are the total 
energies of the incoming (outgoing) particles. The nonrelativistic forms of the cross sections are 
restored when all total energies are replaced by the corresponding masses. 

IV. RESULTS 

In the following subsections the distribution over the breakup phase-space of relativistic effects 
for the exclusive breakup cross section will be presented and comparison to existing data made. 
The five-fold breakup cross sections can be written as 

*» = <„$,' "*""'* "'''*"■ "" 

with the kinematical factor p^in containing the phase-space factor and the initial flux. The dif- 
ferences between the relativistic and nonrelativistic cross sections can result from the dynamical 
part, given by the transition probability for breakup, J2m mr I < 4'o\Uo\4' > Pi or/and from the 
kinematical factor. As a measure of the relativistic effect in a particular complete conflguration of 
the outgoing nucleons we take the quantity 

^Q- \rel ( (S'a \nrel 



( 



A = A(0i, 6*2, (1)12, S) = -'i^^i'^^'aa^^ ^^^ nr^i ^ ^^^'^'>- (^7) 



It should be stressed that relativistic and nonrelativistic kinematics lead to different S-curves in 
a plane of kinetic energies E^^'' — Elf^. This makes the definition of relativistic effects for breakup 
ambiguous and dependent on the procedure which is applied to associate points on the relativistic 
and nonrelativistic S-curves. Those two S-curves can differ significantly in some regions of the 
breakup phase-space and/or at higher incoming neutron energies. Thus a reasonable projecting 
procedure is required, especially when the cross sections change drastically along S-curve. In 
such a case an application of an improper projection method can shift those structures inducing 

10 



artificially values of A. In order to avoid such situations we applied the following procedure to 
associate a point on the nonrelativistic S-curve with a given point on the relativistic one. At small 
values the relativistic and nonrelativistic kinetic energies approach each other and in such a case 
we projected in E^i^ — Elf^ plane along direction perpendicular to the axis El"''' of the detected 
nucleon with smaller kinetic energy. In other regions we projected along direction perpendicular 
to the relativistic S-curve. Such a procedure allowed us in most cases to associate properly points 
on the nonrelativistic S-curve to the points on the relativistic S-curve. 

A. Phase-space distribution of relativistic effects 

In order to study the distribution of relativistic effects over the breakup phase-space we per- 
formed the following investigation for the d(n,nn)p breakup reaction. At given energy of the 
incoming neutron we scanned the whole breakup phase-space and associated all regions with A 
values. In order to provide results which could be of interest for future experiments we restricted 
ourselves to regions of phase-space with cross sections sufficiently large (for given 9i, 62 and (/)i2 
cross sections smaller than 20% of maximal value have been rejected) and kinetic energies of de- 
tected neutrons > 5 MeV. To see how relativistic effects depend on energy we carried through 
such a search at two incoming neutron energies S^°^ = 65 MeV and 200 MeV. In order to locate 
uniquely the phase-space regions where relativistic effects enhance or diminish the nonrelativistic 
breakup cross section we show at each energy two sets of three two-dimensional plots for positive 
and negative sign of A, respectively. The first one is the 9i — 62 plane for the two angles of the 
neutron detectors. The second one is the 9i — <j)i2 plane, where 0i2 = \(l)i — </'2| is the relative 
azimuthal angle for the two detectors. Finally, the third is the Ei — E2 plane for the correlated 
lab. kinetic energies of the two detected neutrons. To fill those three planes we proceed as follow. 
The whole phase-space is filled with discrete points corresponding to certain grids in 9i, 92, 0i, (/>2, 
and El. For 9i and ^2 fixed we search for the maximal magnitude of A (at given sign of A) in the 
three-dimensional subspace spanned by (pi, (f)2, and Ei . Then we combine those maximal A values 
into four groups and associate certain gray tones (colors) to those group values. Next we choose 
a fixed 9i and 0i2 = |<A2| (by putting 0i = 0) and search again for the maximal values of A in 
the two-dimensional subspace spanned by 92 and Ei. The same gray tones and groupings are then 
applied. Finally, in the Ei — E2 plane we search for the maximal A values in the three-dimensional 
subspace spanned by 9i, 92, 4>i2 and repeat the procedure. 

Results of applying that procedure are shown in the ffi'st row of Figs. ^Elfor E^ = 65 MeV 
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and of Figs. 151111 for 200 MeV, separately for positive (Figs. ^ and ^ and negative (Figs. Hj and 
I3J) values of A. Distinguishing between positive and negative A's allows us to locate regions of 
breakup phase-space where relativity increases or diminishes the nonrelativistic cross section. Our 
numbers are based on the CD Bonn potential. Since we use only gray tones we split the variation 
of the quantity A into four groups. Based on the meaning of the gray tones and using the first 
row of Figs. 11141 one can proceed as follows. Choosing a region in the Oi — 02 plane with a black 
tone we know that in the 9i — (j)i2 plane there must exist also black region for the same 9i. This 
allows to read off a certain value of (j)i2- Then the angular positions of the two detectors are fixed, 
which defines the S-curve in the Ei — E2 plane. Along such an S-curve there must be again a black 
region, where one can read off the corresponding range of energies. Choosing for instance another 
combination of tones, like a black one in the 9i — 62 plane, white one in the 9i — (j)i2 plane one 
knows that the S-curve in the Ei — E2 plane lies in the white and maybe gray regions. This should 
explain the use of first row from Figs. HEJ It is seen that relativity can act in both directions, 
increasing or decreasing the nonrelativistic cross section. The magnitude of relativistic effects 
increases with energy of the incoming neutron. Whereas at 65 MeV they approach ~ 14% when 
relativity increases the nonrelativistic cross section and ~ 25% in the opposite case, at 200 MeV 
the corresponding numbers are ss 55% and ^ 60%, respectively. The effects are distributed over 
the entire phase-space. It is seen in Figs. 1 1141 that at both energies (however, at 200 MeV more 
clearly) the large relativistic effects have a tendency to localize in phase-space regions with small 
value of the undetected proton energy and when momenta of the two neutrons are coplanar on 
opposite sides of the beam ((/)i2 ~ 180°). Those geometries are around a region of quasi- free- 
scattering (QFS) condition, where the undetected proton is exactly at rest in the lab. system. 
When relativity increases the nonrelativistic cross section (A > 0) large relativistic effects occur 
at positions of detectors 9i and ^2 around 9i + 92 ~ 75°. For A < this region is shifted to 
9i+92 ^ 100° leading to the following pattern of the nonrelativistic cross section variation due to 
relativity. Starting e.g. at fixed 02 from configuration where relativity increases the nonrelativistic 
cross section and increasing 9i we are led to configurations with decreasing A, resulting finally 
in geometries where relativity diminishes the nonrelativistic cross section. In Fig. Ejwe show this 
pattern at 200 MeV for a number of configurations along S-curve for a fixed angle 02 and changing 
9i. Of course the same is true when exchanging 9i and 02- This characteristic pattern can be 
looked for in experimental data. 

In order to get insight into the origin of maximal deviations we show in the second and the 
third rows of Figs. Ill Jl the values of A calculated from the dynamical and kinematical factor of the 
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cross section, respectively. This is done only for configurations from the first row in those figures 
(geometries with maximal changes of nonrelativistic cross section by relativity). It is clearly seen 
that the effect is predominantly due to the dynamical change of the transition matrix element. 
For localized phase-space regions of the large A values mentioned above the nonrelativistic and 
relativistic kinematical factors are comparable. 

B. Comparison to exclusive breakup cross section data 

At higher energies the exclusive breakup cross sections have been measured below the pion 



production threshold at 65 MeV 



m 



M 
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30|, 156 MeV |3l|, and 200 MeV 32]. In 



Figs. IBl lllI we compare our predictions with the data taken at those energies for some complete 
geometries. In order to investigate importance of the boost we show for each geometry the nonrel- 
ativistic cross section (dotted line) together with three relativistic cross sections corresponding to 
different treatment of the boost. The first result is based on our most extensive approach to the 
boosted potential (approximation of Eq. (|12() : solid line), in the second the boost of the potential 
from the 2N cm. to the 3N cm. is totally neglected (approximation of Eq. (|13() : dashed-dotted 
line), and in the third the k-dependence of the first order relativistic boost correction is omitted 
(approximation of Eq. ((11)) : dashed line) . 

For all presented configurations relativistic effects are clearly visible. Cases of increasing as 
well as diminishing the nonrelativistic cross section by relativity are represented. At 65 MeV (see 
FiglHl) the smallest effects are for configurations b), c), and d) where both nonrelativistic and 
relativistic predictions provide a satisfactory description of data. The largest effects are for two 
QFS configurations e) and f), where the inclusion of relativity reduces the cross section by ~ 9%, 
and for the symmetrical space star (SSS) configuration a) where this effect amounts to ~ 7%. 
The SSS geometry corresponds to the situation where for the arc-length S ~ 30 MeV all three 
nucleons have equal momenta which in the 3N cm. system lie in the plane perpendicular to the 
beam direction. Despite its smallness that relativistic effect leads to a better description of data 
for those three configurations than the nonrelativistic result. It seems to explain the small and up 

to now puzzling overestimation of the 65 MeV SSS cross section data l2a| by modern nuclear forces 

1 I^ 

and can account for the experimental width of these QFS peaks 24]. For QFS configurations e) 
and f) the boost effect is significant in the region of QFS peak. Neglecting boost decreases the 
relativistic cross section by ~ 8%. However, approximating boost by Ea. (|14() leads practically to 
the same result as for the full boost. 
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At 156 MeV (see Figs jT] and |HI) the relativistic effects are significantly larger than at 65 MeV. 
Here among configurations shown there are cases where relativity increases the nonrelativistic 
cross section by up to ~ 30% and cases with diminished cross section by up to ~ 50%. For those 
configurations there are regions of kinetic energies Ei where relativity improves the description of 
data. One can even argue, that the pattern described in subsection IIV Al is supported by some 
of the presented data. Similarly to the situation at 65 MeV importance of the boost depends on 
configuration and neglecting it changes the relativistic cross section up to ~ 15%. It seems that 
the restriction to the approximation of £0. 1)14(1 also here is sufficient for the boost treatment. 

The largest relativistic effects are seen at 200 MeV (see FigslMTTj). They can lead to changes 
of the nonrelativistic cross section increasing it even by up to ~ 40% or decreasing it by up to 
~ 60%, leading predominantly to a better description of data, which , however , is not satisfactory 
in all configurations. Importance of the boost is similar to that at 156 MeV. Also here the pattern 
described in subsection IIV Al is discernible. The data of ref. ,32,] shown in Figs J9llI] seem to be 
shifted in energy Ei with respect to both nonrelativistic and relativistic predictions. Independent 
measurement providing cross sections with smaller error bars on the slope would be desirable. 

V. SUMMARY AND OUTLOOK 

We numerically solved the 3N Faddeev equation for nd scattering including relativistic features 
such as the relativistic form of the free propagator and the change of the NN potential caused by 
the boost of the 2N subsystem. The calculations have been performed at the neutron lab. kinetic 
energies Ej^'^ = 65 MeV, 156 MeV, and 200 MeV. To construct the 3N momentum space basis we 
used the relative momentum of two nucleons in their 2N cm. subsystem together with momentum 
of the spectator nucleon in the 3N cm. system. Such a choice of momenta is adequate for relativistic 
kinematics and allows to generalize the nonrelativistic approach used to solve the nonrelativistic 3N 
Faddeev equation to the relativistic case in a more or less straightforward manner. That relative 
momentum in the two-nucleon subsystem is a generalisation of the standard nonrelativistic Jacobi 
momentum p. We neglected the Wigner rotations of the nucleons spins when boosting the 2N 
cm. subsystem to the 3N cm. frame. As dynamical input we took the nonrelativistic NN 
potential CD Bonn and generated in the 2N cm. system an exactly on-shell equivalent relativistic 
interaction v using the analytical scale transformation of momenta. We checked that in our energy 
range the boost effects for this potential could be sufficiently well incorporated by restricting the 
exact expression to the leading order terms in a q/uj and v/to expansion. At 65 MeV and 200 MeV 
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we performed search for magnitudes and signs of relativistic effects on the exclusive nd breakup 
cross sections over the relevant parts of the breakup phase-space. We found, that depending on 
the phase-space region relativity can decrease as well as increase the nonrelativistic cross section. 
The magnitude of effects rises with the incoming neutron energy. While at 65 MeV the effects 
are rather moderate (~ 20%) at 200 MeV they can change the nonrelativistic cross section even 
by a factor of ~ 2. At that energy relativity leads to a characteristic pattern of the cross section 
variation with 6i and 62- Comparison to existing data seems to support this finding. At 65 MeV 
the inclusion of relativity can explain some discrepancies found in the past between theory and 
data. 

The selectivity of the nd breakup reaction provides opportunity to study different aspects of the 
3N dynamics. Since higher energies seem to be more favorable to study properties of three-nucleon 
forces, precise higher energy exclusive breakup cross sections should be used as a very valuable tool 
to test stringently the incorporation of relativity. Especially, the configurations around the QFS 
breakup geometry due to their large cross sections and insensitivity to the details of nuclear forces 
seems to be favored for this purpose. 
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FIG. 1: (Color online) The regions of the d(n,nn)p breakup phase-space at incoming neutron energy E^^^ = 
65 MeV projected onto the ^i — 02, Oi — (/>i2, and Ei — E2 plains, carrying certain values of A from Ea. H27() as 
indicated in the boxes. In the first row the largest positive values of A, which measure the largest increasing 
relativistic effect in the five- fold breakup cross section ^^ '^^ ^^ , are shown based on the CD Bonn potential. 
The second and the third row shows the values of A obtained only with the dynamical and kincmatical parts 
of the cross section, respectively, for configurations from the first row. The white area belongs either to 
phase-space regions not allowed kinematicaly or to regions rejected by our cuts on the cross sections or 
energies (see text). 

FIG. 2: (Color online) The same as in Fig. ^ but for the largest negative values of A at -E'f'' = 65 MeV 
shown in the first row (the largest decreasing relativistic effect). 



FIG. 4: (Color online) The same as in Fig.|21but for Elf = 200 MeV. 



FIG. 3: (Color online) The same as in Fig. [T] but for E'^'' = 200 MeV. 
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FIG. 5: The five-fold cross section ^^^ ^^ ^g for d{n,nin2)p breakup reaction at El^'' = 200 MeV at fixed 
02 = 37.5°, (/)i2 = 180° and varying Oi: a) Oi = 27.5°, b) Oi ^ 32.5°, c) Oi = 37.5°, d) 9i = 42.5°, e) 
di = 47.5°, f) 9i — 52.5°, g) 6i — 57.5°, h) 9i — 62.5°. The nonrelativistic and relativistic cross sections are 
shown by dashed and solid lines, respectively. The calculations are based on the CD Bonn potential (see 
text). 
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FIG. 6: The five- fold cross section 



for d{n, ni7i2)p breakup reaction at £'„ = 65 MeV for the 



folfowing positions of outgoing neutron detectors: a) 9i = 62 = 54°, (/)i2 = 120°, b) 61 = 20°, 82 — 45°, 
012 = 180°, c) Oi ^ 20°, 92 = 75.6°, </>i2 = 180°, d) 0i = 20°, 02 = 116.2°, 012 = 0°, e) 0i = 30°, 6*2 = 59.5°, 
012 — 180°, f) 6*1 = ^2 = 44°, 012 ~ 180°. The nonrelativistic and relativistic (fuU boost) cross sections 
are shown by dotted and solid lines, respectively. Neglecting the boost totally leads to dashed-dotted line. 
When boost is approximated by Eg. 1141) results in dashed line. All calculations are based on the CD Bonn 

n n n 

potential (see text). The d{p,pp)n data are from: |2y| a) ; |23| b), c) and d) ; [2,91 e) and f). 
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FIG. 7: The five-fold cross section 



d^a 



cin dn dE ^"-"^ '^("'' "■"■)P breakup reaction at £'„ = 156 MeV at fixed 
position 6*1 = 45°, 0i2 = 180°, and varying O^: a) 6I2 = 25°, b) 6*3 = 32.5°, c) 6*3 = 37.5°, d) 62 = 40°, e) 
6*3 = 42.5°, f) 6*3 == 45°, r) 6I2 = 47.5°, h) 62 = 52.5°, i) 6*3 = 60°, j) 62 = 67.5°. Curves as in FigEI The 
d{p,pp)n data are from }31i] . 



20 



> 



C/3 






(N 



a 
a 




30 60 90 30 60 90 30 60 90 

E^ [MeV] E^ [MeV] E^ [MeV] 

FIG. 8: The same as in Fig[7|but at fixed position Oi = 52.5°, </>i2 = 180°, and varying 6*2: a) 62 = 25°, b) 
6*2 = 30°, c) 92 = 32.5°, d) 62 = 35°, e) 62 = 37.5°, f) 6*2 = 40°, g) 6*2 = 45°, h) 62 = 52.5°, i) 6*2 = 60°. 
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FIG. 9: The five-fold cross section 



for d{n, np)n breakup reaction at £'„ — 200 MeV at fixed 



dSlidSl2d£i 

position Oi = 35°, 0i2 = 180°, and varying 6*2: a) 6*2 = 35°, b) 62 = 45°, c) 6*2 = 55°. Curves as in FigEI 
Tlie d{p, pn)p data are from ^32] . 
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FIG. 10: The same as in FigElbut at fixed position 6i = 45°, 0i2 = 180°, and varying 6*2: a) 62 = 35°, b) 
92 = 45°, c) 92 = 55°. 
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FIG. 11: The same as in FigElbut at fixed position Oi = 52°, 0i2 = 180°, and varying 6*2: a) 62 = 35°, b) 
6*2 = 45°. 
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